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Abstract. We study the effects of the confining conditions on the occurrence of stochastic resonance (SR) 
in continuous bistable systems. We model such systems by means of double- well potentials that diverge like 
\x\ q for \x\ — > co. For super-harmonic (hard) potentials with q > 2 the SR peak sharpens with increasing 



q, whereas for sub-harmonic (soft) potentials, q < 2, it gets suppressed. 
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1 Introduction 



The simplest dynamical system displaying stochastic res- 
onance (SR) is a Brownian particle bound into a one- 
dimensional double well under the action of a time os- 
cillating tilt and subjected to fluctuating forces (noise) [U 
12] • The SR mechanism can be revealed as a maximum 
in the amplitude of the periodic component of the aver- 
age particle position as a function of the noise intensity 
(temperature). Due to fluctuations, the particle randomly 
jumps between the two potential wells with Kramers rate 
[3] that depends on the double well potential and temper- 
ature. When the average escape time of the particle out 
of the potential minima (i.e., the inverse of the Kramers 
rate) approximately equals the half time-period of the ap- 
plied perturbation, the noise induced interwell jumps and 
the periodic force synchronize, thus leading to SR. 

When studying the problem of a Brownian particle in 
a symmetric double well periodically tilted in time, the 
corresponding potential U(x) is usually assumed to di- 



2 Model 

The model discussed in the following represents an over- 
damped Brownian particle with coordinate x. Rs dynam- 
ics is described by the Langevin equation, 



r 1 x = -U'{x)+A(t)+£{t). 



(1) 



where (. . .)' = d(. . .)/dx. The confining potential, 

U{x) = U exp(-x 2 /L 2 a ) +k\x\y q , (2) 

is obtained by superimposing a Gaussian repulsive barrier 
of height Uq and width Lq, to a power-law potential well. 
To ensure confinement, our analysis is restricted to q > 1. 
The total potential is mirror symmetric at x = 0, i.e. 
U(x) = U(—x). Depending on q a potential U(x) is called 
hard (super-harmonic) for q > 2, or soft (sub-harmonic) 
for q < 2 [4]. The periodic drive A(t) is chosen as 



A{t) = A cos{f2t) 



(3) 



verge like U(x) 



at large x [213], so as to ensure a 



robust confining action. However, the divergence of the 
potential for |x| — > oo strongly affects the response of the 
system to an external time-periodic forcing. The goal of 
the present paper is to investigate how the Brownian mo- 
tion in a double well changes with the confining strength 
of the one-dimensional potential U(x). For simplicity we 
assume that U(x) ~ \x\ q for x — > ±oo. By studying the 
dependence of a SR spectral quantifier on q, we conclude 
that bistability is a necessary, but not sufficient condition 
for a one-dimensional system to exhibit SR. 



with amplitude, Aq, angular frequency, S7 = 2-7W, and 
time origin arbitrarily set to zero. The fluctuating force 
£(i) is modeled as a stationary zero- mean Gaussian noise 
with auto-correlation function = 2nkBTS(t—t'). 

Here T is the temperature and rj the friction coefficient. 

For numerical purposes it is convenient to choose Uq, 
Lq, and r = i]Lq/Uq as the new units respectively of en- 
ergy, space and time. Correspondingly, the variables and 
the parameters appearing in Eq. |T]) can be replaced by the 
dimensionless quantities x = x/L , t — t/r, k = L^k/Uo, 
A = AqLo/Uq, h = Qr, and f = k B T/U Q . To avoid 
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Fig. 1. Rescaled potential (j4| for k = 0.2 and q ranging be- 
tween 1.5 and 8. The barrier height is approximately constant, 
AU ~ 0.66, and the minima ±x m slowly shrink with q from 
x m ~ 1.59 down to x m ~ 1.17. 



a cumbersome notation, in the following we omit all the 
tildes. In dimcnsionless notation the potential §2$ reads, 

U(x) = exp(-x 2 ) + k\x\ q /q, (4) 

and the Langevin equation {T]) can be rewritten as 

x = 2.xcxp(-x 2 ) - k\x\ q /x + A cos{m) + VT£{t) , (5) 

after the Gaussian noise has been further rescaled so 
that (£(t)) = and (£(*)£(*')) = 2<5(t - t'). In the follow- 
ing we study how changing q influences the response of 
the particle to the periodic forcing signal A(t). As a result 
of rescaling, the height, C/o, and the width, Lo, of the po- 
tential barrier, as well as the friction coefficient, 77, have 
been set to one. The remaining tunable parameter k of 
the potential (|4]) will be kept fixed to k — 0.2 throughout 
the present paper. Due to the Gaussian nature of the po- 
tential barrier, the barrier height, AU , and the potential 
minima, ±x m , weakly depend on q (see Fig.[T]); therefore, 
the observed residual SR dependence on q is mostly an 
effect of the varying confining strength of the potential. 

We have simulated the behavior of the system by nu- 
merically integrating the rescaled Langevin equation |5|) 
through a Milshtein algorithm [HE]. Stochastic trajecto- 
ries were simulated for different time lengths i max and time 
steps At, so as to ensure appropriate numerical accuracy 
and transient effects subtraction. Average quantities have 
been obtained as ensemble averages over at least 10 4 tra- 
jectories. 
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Fig. 2. Rescaled amplitude x(T)/Ao, defined in Eq. ©, 
versus T for the potential © with k = 0.2, q = 2. The 
dashed lines represent the intrawell oscillations, Eq. (JSJ, with 
k = l/|2feln(fe/2)| for T -> 0, and k = k in the limit T -> 00. 



with amplitude, x(T), and phase, 4>(T), approximated by 



A (x 2 



2r 



= — — ■ — — • 

T V4r 2 + fi 2 

4>{T) = arctan(f?/2r) . 



(7) 
(8) 



Here r oc exp(— AU/T) is the Kramers rate and (x 2 )o 
the variance of the stationary unperturbed process x(t) 
(Aq = 0), both temperature dependent quantities. The 
amplitude x(T) can be manipulated by tuning the noise 
level. Note that Eqs. (J6j)- (JSj hold in the linear response 
theory limit, only, i.e., for Aox m <C T and Q > r [HIE]. 

According to Eq. ((7j), in the limit T — > the ampli- 
tude x(T) vanishes due to the potential barrier. The rate 
r for the particle to overcome the potential barrier de- 
creases to zero exponentially when lowering the tempera- 
ture, that is r <C fi. The inter well jumps are thus inhibited 
and the particle gets locked in either minima with prob- 
ability 1/2; hence limr— yo{%) = 0. In contrast, for high 
temperatures, T 3> AU, r may grow much larger than £1 
and, consequently, x(T) ~ (x 2 )q/T. For a hard potential 
with q > 2 we show below that (x 2 )o ~ T 2 l q , so that, 
again, Xvo\t~^oo x(T) = 0. The occurrence of these limits 
for T — > and T — > 00 implies the existence of a maximum 
of x(T) for some optimal T ~ AU. This is the so-called 
spectral characterization of SR PQ . 



3 Results 

In the long time regime, after transient effects subsided, 
the response (x(t)} of a particle moving in a symmetric 
bistable potential U{x) under the action of the signal ([3]) 
with small-amplitude, AQX m <C AU, and low-frequency, 
fl <C U"(x m ), results from the interplay of inter- and the 
intrawell dynamics [T] . On ignoring for the time being the 
intrawell dynamics, the system response at low tempera- 
tures is dominated by its harmonic component [TUT ,8,9, 

(x(t-+oo)) = x(T) cos[f?i - <p(T)} , (6) 



3.1 Harmonic confining potentials 

However, even if the approximate results ([6])- (JSj describe 
correctly the occurrence of SR in most bistable systems, 
Figs. [2]and [3] (q > 1) clearly show that for T — > the am- 
plitude x(T) approaches a non-zero limit x(0) > 0. This 
is a characteristic signature of the intrawell dynamics 
IT2"] . Moreover, for (and only for) q = 2 a similar behavior 
occurs also in the opposite limit T — > 00: the curves x{T) 
attain an horizontal asymptote, see Fig. [2] The coexis- 
tence of these two asymptotes, peculiar to q — 2, strongly 
suppresses the SR peak. 
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The nonzero x(T) limits for T — > and T — > oo 
can be explained by noticing that an overdamped Brow- 
nian particle bound to a generic harmonic potential well, 
U(x) = k(x — xq) 2 /2, responds to the signal (J3]) with am- 
plitude 

x = A Q 



n 2 



(9) 



[Note also that its variance in the absence of forcing (A = 
0) is (.t 2 ) = T/k.} 

In the low temperature limit, T — » 0, the particle de- 
scribed by the Langevin equation (JSJ) is locked in either 
the right or left potential well, where it executes additional 
harmonic oscillations around the corresponding minima 
xo = ±x m prtimn 2 , 1 3 j . Such intrawell oscillations should 
not be mistaken for the interwell dynamics described by 
Eq. ([6|) 9j . Their amplitude is well reproduced by Eq. ((9]) 
with k = U"(±x m ) = |2fcln(fc/2)|. 

In the high temperature limit, T — > oo, the fluctuations 
£ (t) may grow so intense that the barrier of the bistable 
potential (j4|) becomes ineffective; the particle is thus ef- 
fectively confined into a parabolic potential with k = k 
and centered at xq — 0. The amplitude of the periodic 
component of the particle response to the external force 
is then described again by Eq. ([9]) but with k = k. 

For small frequencies the rescaled amplitude x/Aq only 
depends on the curvature of the bistable potential at Xo = 
±x m for T -> 0, x/Aq = l/|2fcln(fc/2)|, and at x = for 
T — > oo, x/Aq — l/k. 

The argument above can be easily generalized for any 
value of q at low temperatures, but it becomes untenable 
in the limit T — > oo, where nonlinearity comes into play. 



3.2 Hard confining potentials 

As anticipated above, at high temperatures the presence 
of the central barrier can be ignored. This implies that for 
T -> oo Eq. simplifies to 



A n 



(* 2 >o 



T 



1 J °°dxx 2 exp{-kx q /qT) 
T J Q °°dx exp(-kxi/qT) 



(10) 



In Eq. (flQ|) we made use of the inequality r ^> Q and 
of the approximate expression Po(x) — A/"exp(— kx q /qT) 
for the stationary probability density of the unperturbed 
process (|5]l: Af is an appropriate normalization constant. 
Note that for sufficiently low J?, the condition r>J] can 
be consistent with the approximations in Eq. 0, whereas 
suppressing the potential barrier makes the very definition 
of r meaningless. 

An explicit calculation yields 



x(T) (q\ 2 /ir(i/q) 



(I) 



1 



A \kJ r(l/q) T 1 - 2 /?' 
Ignoring the algebraic factors we conclude that 



lim x(T) ~ T 



2/9-1 



(11) 



(12) 



From here one can see that x decreases with increasing 
T only for hard confining potentials with q > 2. In par- 
ticular, for the prototypical case of a quartic potential, 
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Fig. 3. Rescaled amplitude x(T)/Ao, defined by Eq. versus 
T for the potential (J4j) with k = 0.2 and different q > 2 (hard 
potentials). The dashed lines are the decay power law T 2//<J_1 . 



q = 4 rjj, one finds x(T) ~ 1/VT, as confirmed by the 
simulation results (see Fig. [3j. For q — 2, one recovers the 
harmonic limit discussed in the foregoing subsection. 

The decay law of x(T), Eq. ([12]) , is clearly a conse- 
quence of the nonlinearity of the potential. Indeed, the 
same power law can be recovered by implementing the 
stochastic linearization scheme of Ref. 14 : In Gaussian 
approximation, for q an integer, ]ha\ x \^ >00 U(x) — kx 2 /2 

with k — (q — l)\\k{x 2 ) q J 2 ~ l ; from the relation (x 2 )q = 
T/k, holding for harmonic potentials, Eq. (TT2]) follows. 

Moreover, x(T) cannot decrease faster than T _1 , which 
happens for q — > oo. It should be noticed that x(T) ~ T _1 
is the decay law predicted in two-state model approxima- 
tion [7], where (x 2 )q is replaced by x^ (i.e., a constant). 



3.3 Soft confining potentials 

Equation (TT2"|) for q < 2 suggests that x(T) may diverge at 
high temperatures. However, when dealing with soft po- 
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Fig. 4. Rescaled amplitude x(T)/Aq versus T for the potential 
((4| with k — 0.2 and q — 1.5 (soft potential). The dashed lines 
represent the horizontal asymptotes 1/1? (see text). In place of 
the SR peak an inflexion point is detectable for low 1? = 2nv. 



is super-harmonic. Most notably, for harmonic and sub- 
harmonic potentials the periodic component of the system 
response may increase monotonically with the noise level. 

(ii) In many experimental reports (see, for a review, 
Ref. [IB]), the authors tried to characterize the SR peak 
by means of Eq. ([6j , without paying much attention to the 
T dependence of the quantity (x 2 )q. In some cases they 
adopted an outright two-state model with (x 2 )o = x^. 
This led to a poor fit of the decaying tail of x(T), whereas 
a more accurate fit could have given a valuable clue to 
better model the system at hand [T?] . 



Acknowledgments 

This work has been supported by the Estonian Science 
Foundation through Grant No. 7466 (M.P., E.H.), Spanish 
MEC and FEDER through project FISICOS (FIS2007- 
60327), and ESF STOCHDYN project (E.H.). 



tentials, the linear theory approximations ([6])- (JS]) must be 
used with caution. In the limit T — > the interwell oscil- 
lation amplitude ([7]) is known to apply only for very small 
perturbation amplitudes [4] : This explains the residual Ao 
dependence of the low T plateaus reported in Fig. |U 

More importantly, in the high T limit, although the 
barrier of a soft potential is awash with noise, confinement 
gets so weak that the particle is driven up and down the 
potential walls primarily by the deterministic force A(t), 
rather than by the noise. [For a comparison, we remind 
that a particle falls from ±oo down to ±x m in a finite time 
for q > 2 and in an infinite time for q < 2.] In conclusion, 
on assuming that the Brownian particle oscillates as if it 
were (almost) free, its amplitude would read 

lim x{T) ~ Aq/Q. (13) 

T — >oo 

x(T) is then expected to develop high T plateaus also 
for q < 2, but, in contrast with the cases discussed in 
Sec. 13. 1\ such plateaus are inverse proportional to the drive 
frequency (also for low frequency drives, see Fig. [4]). 

In the case of sub-harmonic bistable potentials the 
hallmark of SR is thus the monotonic increase of the re- 
sponse amplitude with T, as opposed to the occurrence of 
a maximum often detected in the super-harmonic poten- 
tials. Such a behavior resembles the phenomenon of "SR 
without tuning" discussed in Ref. [15] , with the important 
difference that here it has been observed in a single unit, 
rather than in a summing network of N excitable units. 

4 Conclusions 

We conclude this note with two important remarks: 

(i) The coexistence of two locally stable minima sepa- 
rated by a potential barrier is commonly advocated to ex- 
plain the occurrence of a SR peak in a continuous bistable 
dynamics. Here we have shown that this keeps being true 
as long as the confining action exerted by the potential 
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